Bose-Einstein Condensation in liquid ^He near the liquid-solid transition line 
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We present precision neutron scattering measurements of the Bose-Einstein condensate fraction, 
no(r), and the atomic momentum distribution, n*(k), of liquid ''He at pressure p =24 bar. Both 
the temperature dependence of no(T) and of the width of w*(k) are determined. The no(T) can be 
represented by no(T) = no(0)[l-(T/TA)^] with a small no(0) = 2.80±0.20 % and large 7 = 13±2 for 
T < Tx indicating strong interaction. The onset of BEC is accompanied by a significant narrowing 
of the n*(k). The narrowing accounts for 65 % of the drop in kinetic energy below T\ and reveals an 
important coupling between BEC and fc > states. The experimental results are well reproduced 
by Path Integral Monte Carlo calculations. 



PACS numbers: 03.75.Kk, 78.70.Nx,67.80.bd 

Bose-Einstein condensation (BEC) is pervasive in con- 
densed matter and the origin of spectacular properties 
BEC may be defined as the condensation of a 
macroscopic fraction of Bosons into one single particle 
stated, 13' the onset of long range order in the one- 
body density matrixQ or in a pair function. The phase 
of the macroscopically occupied single particle state or 
pair function introduces phase coherence in the system 
which is the origin of superfluidity and superconductiv- 
ity. Magnetic order is also regularly described^ in terms 
of condensation. BEC in a gas of photons has been 
observed @. Particularly, remarkable properties in dilute 
gases in traps arise from BEC and superfiow. In gases 
the fraction, no, of Bosons in the condensate can be 100 
% and BEC is easier to observe than superflow. In con- 
trast, in dense systems such as liquid "^He, where hq is 
small, superfiow was observed long before BEC 7]. To 
date, BEC is uniquely observed in^helium in the dynamic 
structure factor using neutrons f 



HI 



Reports of possible superfiow in solid helium 12 - 
have stimulated renewed interest in BEC in dense Bose 
systems. Observation of BEC in solid helium would be 
an unambiguous verification of superflow but, as yet, has 
not been observed [l7l - [l9| . To bettei^understand BEC 
in dense systems we have measured |2C| the condensate 
fraction in liquid ^He at low temperature as a function 
of pressure up to solidification, p — 25.3 bar. The full 
atomic momentum distribution, n(k), and especially the 
impact of BEC on n(k) in dense systems, is also of great 
interest. 

In this Letter we report precision measurements of the 
temperature dependence of 7T,(k) and no of liquid ''He 
under pressure p = 24 bar. The measurements were 
made on the ARCS instrument at the Spallation Neutron 
Source (SNS), Oak Ridge National Laboratory (ORNL). 
Path integral Monte Carlo (PIMC) calculations are also 



reported. From the observed n(k) we obtain a Bose- 
Einstein condensate fraction no(T) = no(0)[l — (T/Ta)^] 
with no(0) = 2.80 ± 0.20% and 7 = 13 ± 2 below the 
normal-superfiuid transition temperature Tx— 1.86 K. 
The small value of no and the large value of 7 signal 
strong interaction in the liquid at 24 bar. In addition 
to BEC, the momentum distribution of the atoms above 
the condensate, denoted n*(k), narrows below Tx- With 
the improved precision on ARCS, we are able to deter- 
mine both the temperature dependence of no(T) and the 
width of n* (k) simultaneously. The temperature depen- 
dence of the width of n*(k) below Tx tracks no(T). This 
signals a coupling between BEC and the occupation of 
the higher momentum states. Below Tx, there is both 
BEC and a re-distribution of occupation of the fc > 
states in an interacting Bose system. 

At temperatures close to, but above Tx, the kinetic en- 
ergy, (K), is dominated by quantum zero-point motion 
and hence is relatively temperature independent. When 
the liquid is cooled below Tx, the kinetic energy, (K), 
drops precipitously both due to the onset of BEC and as 
a result of a narrowing of n* (k) with temperature. The 
present results show that at 24 bar approximately 65% 
of the observed drop in (K) comes from the decrease in 
the width of n'^(k) while 35% arises from the onset of 
BEC. This means that determinations of the condensate 
fraction, no, from the drop in (K) below Tx must take 
account of this narrowing of n* (k) . Otherwise no will be 
overestimated. The impact of the narrowing is relatively 
smaller at saturated vapor pressure (SVP) where no is 
larger. However, this effect explains why determinations 
of no from the (K) made assuming a change in weight, 
but no change in the shape, of n*'(k) with temperature 
yield large values of no [10|, ]2m ■ The temperature depen- 
dence of both no{T) and the width of n*(k) are well 
reproduced in PIMC calculations. 
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FIG. 1: Observed scattering intensity S{Q,ijj) as a function 
of energy transfer E — fko and momentum transfer hQ from 
liquid *He at p = 24: bar and T = 40 mK. Signal from the 
empty Al container has been subtracted. The dashed line is 
the calculated *He recoil line, Er = h'^Q^/2m, shown as a 
guide to the eye. 



The atomic momentum distribution is observed in the 
dynamic structure factor, S{Q,uj), at high momentum, 
HQ, and energy, hui, transfer. In this limit, denoted the 
impulse approximation (lA), the energy transfer to the 
sample by the scattered neutrons is quadratic in Q, and 
centered around the ^He recoil line Er = h^Q^/2m, as 
shown in FiglTJ In the lA, S{Q,uj) is conveniently ex- 
pressed in terms of the y-scaling variable, y=(a; — a;j,)/wr, 
yielding[ii[lll, 

JiAiy) = VrS{Q, Lo)= [ dkSiy - fcQ)n(k), (1) 



where kq — and Vr — hQ/m. JjAiy) is denoted 

the longitudinal momentum distribution and its Fourier 
transform, J/a(s), given by Jia{s) = J^^ JiA{y)e~^y''ds 
is the one body density matrix (OBDM) for displace- 
ments s = r.Q along Q. At finite Q, the observed J{Q, y) 
is broadened by final state interactions [2^ and the in- 
strument resolution. Accounting for these effects, single 
particles dynamics such as no, and n(k) are directly ob- 
served from J{Q,y). 

The ARCS instrument was set in its high resolution 
mode and a neutron incident energy Ei = 700 meV was 
selected to allow access to wavevectors up to Q = 28 
A~^. Fig. [1] displays the net 2D contour map obtained 
from liquid "^He after background subtraction. Fig. [2] 
shows the observed J{Q, y) of liquid ''He at Q = 24 
at temperatures below and above T\ along with the in- 
strument resolution function. The relative increased in 
intensity at y = as the temperature is lowered below 
Tx is attributed to the onset of BEC. 

To analyze the data, we follow methods tested previ- 
ously in liquid ^He at SVP d, M, Specifically, 
we express J(Q, s) as a product of the ideal JiAis), and 
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FIG. 2: Observed J{Q,y) at Q = 24 A"^ and at the tem- 
peratures indicated. The dashed line is the measured ARCS 
resolution function at Q = 24 A~^. The increased peak height 
at low temperature is attributed to the onset of BEC. 



the final state (FS) function R{Q,s) [22,l25|. Jia{s) and 
R{Q, s) are determined separately from fits to data. In 
order to extract a condensate no, we assumed as in pre- 
vious work (22l.[25| a model momentum distribution n(k) 
of the form. 



n(k)=no[(5(k)-|-/(k)]-^Ain*(k) 



(2) 



where no(5(k) is the condensate component, n'^(k) is the 
distribution of atoms above the condensate in the k ^ 
states and no/(k) a coupling between the two. The 
Fourier transform of n*(k), Jf^(s) — n*(s), expanded in 
powers of s up to is. 



n*(s) = cxp 



aes" 



6! 



(3) 



The model n(k) in Eq. [5] has thus four adjustable pa- 
rameters, no, 0.2, Oi4, and ae that can be obtained by 
fits to experimental data. Including resolution effects, 
we were able to determine no(T), the momentum distri- 
bution n*(k) and the final state function R{Q,y). 

To get a microscopic understanding of our data, we 
have carried out Path Integral Monte Carlo (PIMC) cal- 
culations of the momentum distribution liquid "^He at the 
same densities and temperatures covered by the exper- 
iment [26|. PIMC is a microscopic stochastic method 
that is able to generate very accurate results relying only 
on the Hamiltonian of the system. The results here pre- 
sented are obtained using a well tested Aziz potential and 
a modern approach based on a high-order action and the 
worm algorithm for a better sampling of permutations 
23. 

Fig. [3]shows the parameters no and 02 obtained by fits 
to the experimental data at several Q values and at T = 
40 mK. The variation with Q arises from the statistical 
uncertainty of the data. The dashed lines indicate the 
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FIG. 3: The parameters a2 and no obtained by fits to data 
at several Q values and temperature T = 40 mK. 



TABLE I: Temperature dependence of the condensate fraction 
no and n(k) parameters in liquid ^He under pressure p =24 
bars . The A transition is at T =1.86 K. The same parameters 
in liquid *He at SVP [25| are shown for comparison. 



P (bar) 


T (K) 


no (%) 


52 (A-^) 


a4 (A-^) 


a6(A-^) 


24 


0.04 


2.88±0.60 


1.10±0.02 


0.63±0.10 


1.35±0.20 




1.00 


2.96±0.70 


1.11±0.02 


0.62±0.10 


1.34±0.15 




1.30 


2.64±0.75 


1.11±0.02 


0.63±0.10 


1.28±0.15 




1.50 


2.56±0.50 


1.12±0.01 


0.61±0.10 


1.26±0.25 




1.75 


1.72±0.70 


1.15±0.02 


0.61±0.10 


1.27±0.25 


24 


1.95 


0.25±1.65 


1.19±0.02 


0.65±0.20 


0.95±0.35 




3.50 


0.32±1.00 


1.18±0.02 


0.56±0.15 


0.98±0.30 




5.00 


-0.04±1.20 


1.18±0.02 


0.53±0.20 


0.44±0.60 


SVP 


0.50 


7.25±0.75 


0.897±0.02 


0.46±0.05 


0.38±0.04 



corresponding average values. The average condensate 
fraction, no, and the parameters 752, 754, and oIq obtained 
by fits to data are listed in Table HI The dependence of 
no{T) and 752 (T) on temperature is shown in Figs. H] 
and El respectively. From Table U and Fig. lU we see 
that no reaches a maximum value of uq = 2.88 ± 0.60 
% at low temperature. As the temperature is decreased 
below T\ = 1.86 K, no{T) increases rapidly toward 
its maximum value. Essentially, noiT) plateaus to its 
maximum value at temperatures close to T\ = 1.86 K. 
This indicates that strong interaction between the ''He 
atoms limits no at higher pressure with a decrease to 
the lowest temperatures unable to reduce np further. 
While the values of the parameters 04 and ae in Table IJ 
are somewhat higher than those observed previously [20| 
at low temperature, the total low temperature n*(k) is 
the same and the value of no at low T is independent of 
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FIG. 4: Observed condensate fraction as a function of tem- 
perature. The dashed line is a line fit to the experimental 
data using of no(r) = no(0)[l - (T/Ta)^], where Ta = 1.86 
K. The open triangles are the simulated PIMC results. 



which n*(k) is used. 

In a Bose gas, no(T) = no(0)[l - (T/Tx)^] with no = 
100 % and 7 = 3/2. A fit of this expression to the ob- 
served no(r) in hquid ^He at SVP ^ gives no(0) = 7.25 
± 0.75 % and 7 = 5.5 ± 1.0. A fit of the same expression 
to the present observed no(T) at 24 bar gives no(0) = 
2.8 ± 0.20 % and 7 = 13 ± 2.0. The fit is shown as a 
dashed line in Fig. 01 The large value of gamma reflects 
the strong interaction in liquid "^He at 24 bar. 

From Table [J and Fig. O we see that the parameter 
752 — (I^Qp) which sets the width of n'^(k) decreases 
from 1.18 A"^ in the normal phase (T > Ta) to 1.10 
-1.11 A^^ at low temperature. That is, while 752 is 
approximately independent of T in the normal phase, 
752 drops abruptly at temperatures immediately below 
T\. This abrupt decrease is unlikely to be a thermal ef- 
fect since the thermal energy ksTx is already much less 
than the zero point energy (approximately {K) = 21.47 
K). Rather, the abrupt drop of 752 below Tx suggests a 
link to the onset of BEC. To test this picture we show the 
function a2(T) = a2{Tx) - A [1 - {T/Tx^] where a2(TA) 
= 1.18 A-2, A = a2{Tx) - 02(0) = 0.075 A-^ and 7 = 13 
as a dashed line in Fig. [5] which has the same temper- 
ature dependence as no(T). The dashed line reproduces 
the observed a2{T) well. Below Ta, there appears to be 
a coupling between no and n*(k), perhaps of the same 
form as no/(k), which leads to a narrowing of n*(k). 

The narrowing of n*(k) below Ta is reproduced by 
PIMC calculations. The present PIMC values of 0.2 are 
shown in Fig. [5] and they also decrease abruptly below 
Ta . Thus the sharp reduction of 752 below Ta , not ob- 
served previously but observable with the increased pre- 
cision of the ARCS neutron scattering instrument, is sup- 
ported by accurate PIMC calculations. Below Ta there 
is both BEC and a narrowing of n'^(k) in a strongly 
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FIG. 5: Temperature dependence of the width a.2 of n*(k) 
at p =24 bar: simulation and experiment. The dashed Une 
shows that a.2(T) has a temperature dependence that tracks 
no(T). The inset shows the corresponding {K){T). 



interacting Bose liquid. 

In summary, the observed and PIMC values of 
no{T) in liquid ^He at 24 bar near the solidification 
line {p = 25.3 bar) are well represented by no(T) = 
no(0)[l - (T/Ta)^] with no(0) = 2.80 ± 0.20 % and 
7 = 13 ± 2. In a Bose gas 7 = 1.5 and in liquid ''He 
at SVP 7 = 5.5 ± 1.0. The large 7 at 24 bar indicates 
strong interaction in the liquid with hq saturating to a 
small value at temperatures close to T^. On cooling be- 
low Tx, there is both BEC and a narrowing of the atomic 
momentum distribution of the atoms above the conden- 
sate, n*(k). The narrowing is characterized here by a 
drop in the width, 0:2 = (I^qI^), of n*(k). The tempera- 
ture dependence of a2 below Tx tracks r7,o(T) indicating 
that interaction between the condensate and the higher 
momentum states causes the narrowing of n*(k). This 
coupling between no and n*(k) is currently not under- 
stood. 

Both BEC and the narrowing of n*(k) contribute to 

the drop in the {K) at temperatures below Tx- The 
present observed and PIMC values of the temperature 
dependence of no{T) and (K) at 24 bar agree well. At p 
— 24 bar, approximately 65% of the drop in (K) arises 
from the narrowing of n*(k) below Tx- Thus an no ob- 
tained from the {K) assuming no narrowing of n*(k) 
would significantly overestimate no- Indeed, if we ap- 
ply the method to the present data, we get no=7.5% at 
40 mK, which is more than 2 times the observed value. 
In the liquid at SVP where rio is larger, the relative re- 
duction of the {K) arising from the narrowing of n* (k) is 
smaller. However, it is still significant and existing val- 
ues of no at SVP determined from the (K) assuming no 
narrowing of n*(k) may have to be corrected. 
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